Quasi dynamical symmetry in an interacting boson model phase transition 
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The oft-observed persistence of symmetry properties in the face of strong symmetry-breaking 
interactions is examined in the SO(5)-invariant interacting boson model. This model exhibits a 
transition between two phases associated with U(5) and 0(6) symmetries, respectively, as the value 
of a control parameter a progresses from to f . The remarkable fact is that, for intermediate values 
of a the model states exhibit the characteristics of its closest symmetry limit for all but a relatively 
narrow transition region that becomes progressively narrower as the particle number of the model 
increases. This phenomenon is explained in terms of quasi-dynamical symmetry. 



There have been numerous recent studies of phase 
transitions in nuclear models 0, 0, HI 01 • Being fi- 
nite zero-temperature many-particle systems, nuclei do 
not exhibit the sharp phase transitions observed in con- 
densed matter physics. Nevertheless, theoretical models 
designed to describe nuclei for particular values of their 
parameters can be extended to study their behavior as 
their parameters arc varied, e.g. as the particle number 
goes to infinity. 

This letter focuses on a phase transition in the in- 
teracting boson model (IBM) of Arima and Iachello 0. 
The IBM comprises N boson particles each of which has 
two states: an L = (s-boson) state and an L = 2 
(d-boson) state with five orientations. The creation 
and annihilation operators for these bosons are denoted 
{s\o%;v = 0,±I ± 2} and {s,d v ;u = 0,±I ± 2}, re- 
spectively. They satisfy the usual boson commutation 
relations 

= i, [^<£] = W*.4] = [dV] = o. (i) 

Thus, the Hilbert space of the IBM carries an irreducible 
representation (irrep) of the group U(6) and can be re- 
alized as a subspace of states of N quanta of a six- 
dimensional harmonic oscillator. The group U(6) has 
several subgroups and different phases of the model can 
be associated with Hamiltonians that arc invariant under 
different subgroups. We consider the Hamiltonian 



H(a) = (1- a) Hi + aH 2 



(2) 



with control parameter a, where Hi is the U(5)-invariant 
d-boson number operator and Hi is the 0(6)-invariant 
operator: 

Hi=n = Y,dU\ H 2 = ^S+S-, (3) 



where 



S + = i(dt-dt_ s t s t) j §_ = i(d-d-na). (4) 

The electric quadrupole moment operator is represented 
in the model by 

Q„ = -=(dts + sU v ) , i/ = 0,±l,±2, (5) 
v ./V 



where Z is a suitable norm factor (can be thought of as 
the charge). 

This letter shows the behavior of the energy-level spec- 
trum and the electric quadrupole transitions for the 
above model as a is varied over the range < a < 1. 
The model is analytically solvable in its U(5) (a = 0) 
and 0(6) (a = 1) symmetry limits. The numerically 
computed results given below for < a < 1 are deter- 
mined by use of the simple SU(1, l) s x SU(1, l) d spectrum 
generating algebra with basis elements 



S — 2" 5 S , 



= id* -dK 



\d ■ d. 



(6) 
(7) 



Energy levels, labelled by an SO(5) quantum number v, 
and E2 transition rates for decay of the first excited state 
of the model are shown as a function of a in figs. 1 1121 




FIG. 1: Excitation energies of the Hamiltonian H(a) plotted 
as a function of a. Precise, numerically computed, energies 
are shown as continuous lines. The dotted lines are the results 
of an RPA calculation (for a < 0.5) and of another harmonic 
aproximation (for a > 0.5). 

A superficial look at the results of figs. ^ an d0 would 
seem to suggest that the model holds onto its U(5) sym- 
metry for < a < 0.3 and to its 0(6) symmetry for 
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FIG. 2: B(E2) transition rates for decay of the first excited 
v = 1 energy level to the ground state in the IBM for various 
values of N expressed in units such that B(E2; 1 — > 0) = 100 
in the U(5) (a = 0) limit. The continuous lines for a < 0.5 
are for the RPA and those for a > are for the harmonic 
approximation. 



0.8 > a < 1 and that there is a transition between the 
states of one symmetry to the other in the intermediate 
region. It turns out this is an overly simplistic view of 
what happens. Insight into the actual evolution of the 
model states with increasing a is given by approximate 
solutions which predict a phase transition and do so with 
increasing accuracy as N increases. We consider the fa- 
miliar Random Phase Approximation for < a < 0.5 
and another harmonic approximation for 0.5 < a < 1. 

For small values of a, the RPA gives quasi-boson exci- 
tation operators 



Dl = (xdts-ysU v )/VN, 



(8) 



with coefficients that satisfy an eigenvector equation of 
the non-Hermitean form 



A B 

-B -A 



= e 



(9) 



with submatrices given, e.g. in the double-commutator 
equations of motion formalism Q , by 



= (1 - (10) 



= -\a8^ v , (11) 



'fiV — N 

where \<f>) is the uncorrelated ground state given by the s- 
boson condensate \<f>) = (s^) n \0)/Vn\. The RPA energy 



spectrum is shown in fig.njas dotted lines for a < 0.5. It 
predicts a collapse of the excitation energies to zero and, 
in fig. El a divergence of the E2 transition rate for decay 
of the v = 1 first excited state as a — ► 0.5. Thus, the 
RPA predicts a phase transition at a cr ;t = 0.5 in accord 
with Thouless' Hartree-Fock stability condition ||. 

Details of the calculations will be given in an arti- 
cle to follow. It will then also be shown that there is 
another harmonic approximation complementary to the 
RPA which gives the spectrum accurately in the a < 1 
region. The essential ingredients of this approximation 
are shown in fig. [3] which gives the expansion coefficients 




FIG. 3: Wave functions for the N — 60 ground state and first 
excited v = state for a = 1 and 0.75. The dots give the 
coefficients of the states in a basis of d-boson number n. The 
continuous lines are those of the harmonic approximation. 

of the lowest and first excited v = states for a = 1.0 
and 0.75 (cf. states labelled by v in fig. ^) in a basis 
{|nu=:0)} which diagonalizcs the Hamiltonian Hi (n is 
the (i-boson number). The smooth curves through these 
numerically computed coefficients are simply harmonic 
oscillator wave functions centred about N/2 for a = 1 
and shifted to smaller mean values of n, in the manner 
of coherent state wave functions, as a is decreased. The 
energy levels and E2 transition rates predicted by this 
approximation are shown for 0.5 < a < 1 in figs. 1 1121 
They are precise for a = 1 and become increasingly ac- 
curate for smaller a with increasing values of the boson 
number N. 

An examination of the physical significance of the wave 
functions of fig. |21 reveals that they represent a nucleus 
with a large mean quadrupole deformation at a = 1 that 
decreases as a falls and the nucleus moves towards a 
spherical equilibrium shape. Such behavior is expected 
from the form of the classical potential energy corre- 
sponding to the Hamiltonian @. This potential, given 
as a function 

V Q (/3)=iV[/3 2 +a(-2/3 2 +/3 4 ) + 0(l)] (12) 

of a and a classical quadrupole shape variable /3, is shown 
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FIG. 4: The classical potential energy V a (/3) for different val- 
ues of a. 

for a few values of a in fig. The magnitude of the E2 
(electric quadrupolc) gamma-ray transition rate shown in 
fig-Elis closely related to the mean square quadrupolc mo- 
ment of the ground state and a good indicator of how the 
mean value of the quadrupolc moment actually evolves 
with a in the model. 

Figure [3] reveals that the eigenfunctions of H(a) for 
0.75 < a < 1 are close approximations to harmonic oscil- 
lator (Glauber) coherent state wave functions. However, 
as a falls below 0.75 (for N = 60), the wave functions 
reach the lower n = boundary and the harmonic ap- 
proximation begins to break down. For larger values of 
N, when the potential V a becomes deeper and the widths 
of the harmonic oscillator wave functions narrower, the 
approximation holds for smaller values of a. However, 
at a = 0.5, the centroid of the wave functions of the 
harmonic approximation lie precisely at n = 0; thus the 
model breaks down asm 0.5 for all values of N. This 
is what one would expect for the quantum states of a 
model with potential energy given by eqn. (|12|l . 

Similar insights can be gained for small a from the 
RPA. For large values of N, in which the RPA is equiv- 
alent to the Bogolyubov approximation of replacing the 
s-boson operators and s by the c-number y/N, the 
RPA transformation of eqn. iJSJ is seen as an SU(1,1) 
transformation 

4^ Dl = xdt-yd v . (13) 

The corresponding transformation of the ground state (to 
a D-boson vacuum) is then a transformation to a dilated 
(anti-squcczed) coherent state in which the equilibrium 
quadrupole shape remains spherical but the width of its 
(harmonic oscillator Gaussian) wave function, along with 
the magnitude of the quadrupole shape fluctuations grow 
with increasing a. As a approaches the critical value 0.5 
from below, the spherical equibrium shape, indicated by 
the potential shown in fig. becomes unstable and the 
width of the D-boson vacuum wave function diverges. 

The precision of the RPA and the other harmonic ap- 
proximation, for low-energy model states outside of a 
transition region that becomes increasingly narrow as 
N — > oo, suggests a useful definition of the concept of 
phase in such situations. 



Definition (phase): If SGC is a subgroup chain 

Gx D G 2 C • • • (14) 

of transformations of a system, then a subset of states of 
the system is said to be in an SGC-phase if the proper- 
ties of the subset are indistinguishable (to within a spec- 
ified accuracy) from a corresponding subset of states of a 
model whose eigenstates belong to irreps of the subgroups 
in the chain and whose obscrvablcs arc infinitesimal gen- 
erators of G\. 

The above results show that there is domain of a ex- 
tending from to some upper limit below 0.5 for which 
a subset of low-energy states are described accurately by 
the RPA. The RPA formalism shows that these states are 
cquivalently described by a U(5)-invariant model Hamil- 
tonian 

Hkpa = e ^2 D i D " + constant , (15) 
with electric quadrupole operator 

Q« PA = ^=(^t S + S t^ ); ^ ,±l,±2. (16) 
V N 

Thus, according to the definition, the subspace of states 
that are described to within the required limits of ac- 
curacy by the RPA, for a given value of a, belong to a 
U(6) D U(5) D SO(5) phase. 

It is similarly shown that the subset of states whose 
energies cluster around those of the harmonic approxi- 
mation for a < 1 span an 0(6) irrep (for some choice 
of 0(6)). Moreover, to the extent that the harmonic ap- 
proximation gives an accurate description of a subset of 
states in this region, the energy levels and E2 transitions 
between these states are accurately described by a model 
whose states are labelled by the quantum numbers of a 
U(6) D 0(6) Z> SO(5) subgroup chain (the small spread 
of energies of a cluster is readily accommodated by in- 
cluding a term proportional to the SO(5) C 0(6) Casimir 
invariant in model Hamiltonian). Thus, by the definition, 
the states that are accurately described by the harmonic 
approximation, belong to a U(6) D 0(6) D SO(5) phase. 

The above analysis of the nature of the solutions of the 
IBM Hamiltonian (J2J in terms of the RPA shows why, in 
spite of the fact that the U(5) symmetry of H (0) = Hi is 
quickly broken by the a7?2 term in H(a) for a > 0, the 
results look as though the U(5) symmetry is preserved. 
Similarly, the harmonic approximation of 0.5 > a < 1 
shows why the apparent 0(6) symmetry is also retained 
over a considerably larger region of a than expected. 

An apparent persistence of a symmetry when numer- 
ical calculations show there to be strong mixing of the 
irreps of the symmetry group in question has been ob- 
served many times. It has been expressed in terms of 
what has been dubbed quasi- dynamical symmetry |l0j |: 
cf. rcf. ;llj for a review. 

A sense of what quasi-dynamical symmetry means is 
obtained by recalling that different but equivalent irreps 
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of a group can be combined coherently to give new (equiv- 
alent) irreps. A simple example would be the coherent 
mixing of the basis states {\aLM)} and {\bLM)} of two 
irreps of SO (3) of angular momentum L to form new basis 
states 

\cLM) = C a \aLM) + C b \bLM) (17) 

for an equivalent mixed irrep of the same L; coherent 
mixing means that the coefficients C a and Cb are in- 
dependent of M . Quasi-dynamical symmetry arises be- 
cause many groups have distinct irreps that are similar 
and difficult to distinguish, especially by consideration of 
a subset of their states. Thus, for example, it is possi- 
ble to represent many states of a free particle with wave 
functions of good linear momentum even though the real 
states are wave packets comprising coherent mixtures of 
plane wave functions of similar momentum. 

It is not surprising then to find that subsets of states 
of a given system can be described by models with dy- 
namical symmetries that are only quasi-dynamical sym- 
metries of the sytem being described. Because a model 
can at best describe a subset of the states of a real phys- 
ical system to within limited accuracy, it could not be 
otherwise. However, an explicit recognition of the pos- 
sible quasi-dynamical symmetries of physical systems is 
invaluable for interpreting the implications of successful 
models and in designing successful approximations. In 
this note, I have explicitly determined realizations of the 
quasi-dynamical symmetries of two phases of the IBM 
and thereby obtained an explanation of why these sym- 
metries appear to persist in spite of the known presence 
of strong symmetry mixing interactions. Simply stated, 
the mixed states of the original dynamical group become 
the unmixed states of a quasi-dynamical group. 



The above perspective leads one to think of the evolu- 
tion of the low-energy states of the IBM as they progress 
though a phase transition in terms of the evolution of the 
quasi-dynamical group. In approximate pictorial terms, 
the effect of the symmetry-breaking interaction is pri- 
marily to distort the quasi-dynamical symmetry, rather 
than break it until a critical point is reached at which it 
can be distorted no more. At this point the symmetry 
really starts to break up; the system enters the transi- 
tion region and, as it emerges on the other side, a new 
quasi-dynamical symmetry associated with another other 
phase begins to develop. 

An interesting characteristic of the above results is that 
for reasonably large values of the boson number N, the 
low-energy states of the model can be assigned to one of 
three domains: one for which the U(6) D U(5) D SO(5) 
phase is appropriate, one for which the U(6) D 0(6) D 
SO (5) phase is appropriate, and an intermediate transi- 
tion domain which shrinks with increasing N. A recent 
suggestion of critical point symmetries which apply in 
the middle of the transition region is therefore of consid- 
erable interest; for example, it suggests that it might be 
meaningful to think of an intermediate phase separating 
the other two. This suggests that an interesting sequel 
to the present investigation would be an exploration of 
the ways the properties of the system in the transition 
domain evolve function of the boson number N . 
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